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Abstract. - Force-induced reversible transformations between coiled and normal polymorphs of 
bacterial flagella have been observed in recent optical-tweezer experiment. We introduce a discrete 
elastic rod model with two competing helical states governed by a fluctuating spin-like variable 
that represents the underlying conformational states of flagellin monomers. Using hybrid Brownian 
dynamics Monte-Carlo simulations, we show that a helix undergoes shape transitions dominated 
by domain wall nucleation and motion in response to externally applied uniaxial tension. A 
scaling argument for the critical force is presented in good agreement with experimental and 
simulation results. Stretching rate-dependent elasticity including a buckling instability are found, 
also consistent with the experiment. 
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i: Introduction — Many motile bacteria achieve directed 
jjropulsion in a viscous environment by rotating flagel- 
lar filaments [1]. A number of common bacteria, such as 
Escherichia coli or Salmonella, are externally flagellated, 
where each flageUum consists of a helically-shaped elastic 
filament that is attached to a rotary motor embedded in 
the rigid cell wall via a flexible hook. The flagellar flla- 
ihent is a large homogeneous assembly of a single protein, 
flagellin, but transforms between several helical shapes of 
different pitch, radius and handedness in response to envi- 
Eonmental or mechanical stimuli such as change of temper- 
ature, pH, or externally applied forces and torques [2]. For 
example, Salmonella swim by rotating their left-handed 
helical filaments counter-clockwise, forming a co-ordinated 
bundle of filaments [3,4]. The reverse of the motor rota- 
tional direction gives rise to a tumbling motion of the cell 
body, often (but not always) accompanied by changes in 
filament chirality to right-handed [2,5]. By regulating the 
relative durations of swimming and tumbling modes, bac- 
teria are able to migrate through a stimulus gradient in 
the favorable direction. 

The flagellar filament is composed of eleven protofila- 
ments that each consist of a stack of flagellin monomers 
and that are wrapped around each other [6] . As flrst pro- 
posed by Asakura [7] and later proved by crystallogra- 
phy [8], a flagellin monomer can take two conformations, 
L-type and R-type, that are slightly different in intrinsic 
twist and length. Each protofilamcnt can be assumed to 



consist entirely of either L or R monomers, however, the 
complete filament is thought to contain variable numbers 
of L and R-type protofilaments. Depending on how many 
protofilaments are of the L or R type, the fiagellum takes 
different polymorphic helical states, thereby minimizing 
the elastic strain energy. A simple geometric model de- 
veloped by Calladine in line with this idea [9] explained 
the observed spectrum of flagellar polymorphic forms. A 
coarse-grained continuum rod model has also been devel- 
oped recently [10]. 

The Asakura-Calladine model has revealed the design 
principle of bacterial flagellar fllament. Recently, dy- 
namics of flagellar polymorphic transitions have moved 
into the focus. Turner et al. have reported an in- 
tricate sequence of transformations of fluorescently la- 
belled E. coli flagella during a single run-to-tumble mode 
switching [11]. Clearly, understanding elastic energy 
barriers between polymorphic forms is crucial for the 
understanding of switching dynamics. Several studies, 
both experimental and theoretical, have addressed this is- 
sue [2,12-14]. Hotani reported cyclic chirality transfor- 
mations of Salmonella flagella tethered at one end to a 
glass surface and subjected to an external fluid flow [12]. 
Above a certain critical flow velocity (a few /zm/s), chi- 
rality flips occur at the tethered end and nucleated chi- 
rality domains grow and propagate steadily down the fll- 
ament. The magnitudes of pulling force or mechanical 
torque necessary for the onset of the normal-to-semi-coiled 
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transition were estimated as 15 pN and 10~^* N m respec- 
tively, based on a slender-body hydrodynamics analysis of 
the experimental geometry. Coombs and Goldstein [14] 
theoretically analysed and rationalized Hotani's observa- 
tion based on a bistable helix model introduced earlier by 
Goldstein and co-workers [15]. Using an optical tweezer 
setup, Darnton and Berg have recently measured the force- 
extension curves of Salmonella flagella filaments [16]. Me- 
chanical stretching induces a reversible polymorphic trans- 
formation between coiled and normal state at a thresh- 
old force around 3-5 pN, characterized by large hystere- 
sis during pulling and relaxing. The measured force is 
shown to be strongly rate-dependent, suggesting that the 
observed polymorphic transformation is kinetically con- 
trolled. Stretching experiments of Salmonella flagella have 
also been performed using atomic force microscopy by 
Cluzel et al. [17]. Since atomistic simulations are impos- 
sible with present computer power due to the huge size 
of a flagcllum and the large time scales involved, coarse- 
graining is imperative. In one recent simulation study, a 
flagellin monomer is represented by 15 point niassc;s with 
interactions parameterized based on the known protein 
crystal structure [18]. Using an implicit solvent model, 
different helical structures arc; observed depending on ro- 
tational direction and in qualitative accord with experi- 
ments. But still, the numerical simulation effort only al- 
lows to simulate the model flagellum for about 30 /is which 
is short compared to intrinsic time scales of polymorphic 
dynamics [18] . It is therefore desirable to take the coarse- 
graining to an even higher level. 

In this paper, we introduce an elastic helix model that 
is coupled to an internal discrete variable that locally de- 
scribes polymorphic states. In the present formulation we 
consider an Ising-like spin variable that switches between 
two states of different helix pitch and radius. We introduce 
a hybrid Brownian-dynamics Monte-Carlo simulation and 
show that our model reproduces well the observed force- 
extension relationship of flagellar filaments accompanying 
polymorphic transformations [16]. As in experiments, the 
force-extension curves are strongly stretching-rate depen- 
dent. We also give simple scaling arguments for the critical 
switching force, in good agreement with experimental and 
our numerical data. 

Model — Within linear elasticity theory, the bending 
and twisting energy for an isotropic helical rod of contour 
length L parameterized by arclength s is given by 

E,i = 1 Tds [ASll + A(02 - S^lf + C{ns - Cl^f] , 
^ Jo 

(1) 

where A and C stand, respectively, for the bending and 
twisting modulus [19]. ft = (Qi, ri2, ^^3) is the strain rate 
vector; the curvature k satisfies = fl\ + fl\ and ^3 
is the twist density. The ground-state shape of a fila- 
ment is specified by the intrinsic curvature kq and tor- 
sion To, related to the radius R and pitch P via kq = 
Att^R/{P^ + 47r2i?2) and tq = 27rP/(p2 + Att^R'^). Al- 



ternatively, one can specify the pitch angle tp through 
t&ntp = 2ttR/P, and the contour length of a single he- 
Ucal turn, £ = (P'^ + Att^R^Y^^. The relation between 
these two descriptions is given by kq = (27r/£)sin^ and 
To = (27r/^) cos V'- Experimentally, the so-called "coiled" 
state of a filament is the equilibrium shape at pH 4 and 
3° C (see the phase diagram in ref. [16]). We in partic- 
ular study the polymorphic transition from the coiled to 
the so-called "normal" state. For the coiled and normal 
forms of Salmonella, which we respectively label as 1 and 
2, £ and ip were measured as (^1,^2) = (3.4 /^m, 2.5 /^m) 
and {tpi,ip2) = (76.7°, 31.3°) [16]. In the present study, 
we for simplicity set the length £ for both forms equal and 
choose pitch angles ipi = 73.3° and tp2 = 29.7°, close to 
the experimental values. The intrinsic curvature and tor- 
sion for each helical state follow as Kom = (27r/£) sin ■!/;„, 
and Tom = (27r/^) cos V'm, where m = 1, 2. 

In the simulations we consider a chain of A^+1 connected 
spheres of diameter a. Each bead is specified by its posi- 
tion Tj and rotation angle about the local tangent , from 
which the discrete strain rate Jlj is calculated. Analogous 
to the helix-coil model for biopolymers [20], each sphere 
carries a discrete variable aj, where cTj = 1 specifies the 
coiled form and aj = — 1 the normal form. To account for 
cooperativity, induced by elastic strain at the interface be- 
tween two polymorphic states, a short-range interaction J 
between neighbouring "spin" variables aj and (Jj+i is in- 
cluded. The total energy functional of the system (apart 
from a stretching contribution introduced later) reads in 
discretized form 



N 



E[fl,a] = — J (TjCJj — ha aj 
(hi) J=i 



(2) 



N 



where («, j) implies summation over nearest neighbors and 
the (T-dependent intrinsic strain rates are fl2j — <^j{i^oi — 
'«02)/2-f (kqi -}-Ko2)/2, and il!]^- = crj(Toi - to2)/2 + (tqi + 
To2)/2. The chemical potential bias per length, h, accounts 
for a preference of one state over the other, and depends on 
environmental conditions such as solution pH or salinity. 

Effective potential — To gain insight into the cou- 
pling between helix elasticity and the polymorphic degrees 
of freedom, we formally define an effective elastic energy 
density f{ft) by integrating out the variables a as 



exp 



N 



-Ein,cr]/kBT 



(3) 



where ksT denotes thermal energy. For uniform helical 
shape, i.e., j-independent strain rates, Cl2j = ^2 and 
flsj = fls, the summation over a can be performed us- 
ing transfer matrix method. Assuming periodic boundary 
condition ai = crjv, and discarding unimportant constants. 
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Fig. 1: Plot of the rescaled effective elastic free energy / = 
fa/ksT, eq. (4), as function of curvature and twist strain rates 
$12, fis (with fli — 0). The two local minima correspond to the 
normal and coiled form as indicated. Parameters values are 
A/{akBT) = 1844, C/A = 0.8, J/fcsT = 10, and h/ksT = 7.7. 

we obtain 

see eq. (4) 

where we defined ksT A{n) = ha + ^{kqi — ko2)[^^2 — 
i(Koi + K02)] + ^(toi -To2)[f^3 - ^(7-01+7-02)] which acts 
as a strain-dependent external field. The effective elastic 
energy /(Jl) is plotted in fig. 1 as a function of the two 
relevant strain rates for one specific parameter set. The 
potential / is a double-well with two minima correspond- 
ing to the normal and the coiled form. Previous theoretic 
treatments used a very similar shape of the elastic free 
energy as a starting point [14]. 

Simulation details — In our hybrid simulations, the 
energy functional in addition includes a stretching contri- 
bution that ensures connectivity of the spheres, E^t = 
~ — a)^. The local elastic transla- 
tional force and torque about the local tangent Tj act- 
ing on each sphere are calculated using the variational 
method described previously [21,22], leading to the cou- 
pled Langevin equations: 

dtVi = fioFi + ^i{t), and dt(j)i = MrTj + ^^(i), (6) 

where dt is the time derivative. Neglecting hydrodynamic 
effects, we take the mobility matrix to be diagonal and 
use the Stokes translational and rotational mobilities of 



a sphere /^q — l/{3Trria) and fj,r — \/{'Krio?), where rj is 
the solvent viscosity. The vectorial random forcings (,{t) 
and 2(t) model the coupling to a heat bath and obey the 
fluctuation-dissipation relations. Standard Monte-Carlo 
Metropolis methods are applied to update a. In each 
Brownian dynamics time step, 20 Monte-Carlo steps are 
carried out, meaning that equilibration of the a degrees of 
freedom is much faster than filament shape relaxations. 

For the simulations we discretize the Langevin equations 
(6) with time step A and rescale all lengths, times and en- 
ergies; the rescaled time step A = AfcsT^o/o^ is for suf- 
ficient numerical accuracy chosen in the range A = 10~^- 
10"^. The stretching modulus is set to K/ksTa^ = 10"* 
which keeps bond length fiuctuations negligibly small. Ob- 
servables are calculated every 10^-10"' steps, total simula- 
tion times are in the order of 10^-10^ steps. One filament 
end is fixed at the origin, and the other end (initially be- 
ing at its equilibrium position, i.e., z = Lcosipi) is moved 
along the z-axis (identical to the helix axis). We set the 
maximum extension z — 0.8L, the total stretching dis- 
tance is thus D — (0.8 — cos'ipi)L. We fix the number 
of beads L/a = N ~ and the arclength of one helical 
turn € = 15a, the pitch angle ■;/; thus uniquely determines 
the helical shape of the filament. We set the bend/twist 
ratio F = C/A = 0.8 with A = A/iakeT) = 1844. The 
bend persistence length is thus Lp/L ~ 31, comparable to 
the experimental conditions Lp/L c± 40 — 120, obtained 
from measured values of the bending stiffness v4 ~ 3.5 
pN-/im^ and typical filament length of 7.6-19.5 /im [16]. 
Although the magnitude of J is experimentally unknown, 
we assume J/ksT 10 sufficiently larger than thermal 
energy in accord with the observation that thermally as- 
sisted nucleation of polymorphic transitions is absent in 
the stress- free state [14,15]. The rescaled stretching speed, 
V = Va/ (fiokBT), is changed in the range V = 0.01-0.2, 
to study systematically the rate-dependent elasticity. As- 
suming a ~ 126 nm (giving L — aN ~ 7.6 fiia for = 60 
as in the experiment [16]), we obtain 0.3 — 3.3 /im/s compa- 
rable to the experimental V ~ 0.4 ^m/s [16]. Our filament 
is thus thicker than a real bacterial flagellar filament by 
roughly a factor of 5, which however is not serious as the 
filament thickness enters hydrodynamic drag coefficients 
only logarithmically. The pulling end of the filament is 
free to rotate, i.e., no external torque is applied, while 
the tethered end is not allowed to rotate, suitable to the 
experiment condition. 

Force-induced shape transformations In fig. 2, 

we show one typical force-extension curve obtained at 
pulling rate V = 0.0124 (the slowest considered), with 
several snapshots of the filament at the points specified in 
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Fig. 2: (a) Typical force-extension curve obtained for N — 60, 
A = 1844, C/A = 0.8, J/fcsT = 10 and ^ = ha/ksT = 7.7 at 
a stretching speed V = 0.0124. Red line is the stretching, the 
green one is the relaxing curve. The broken lines follow from 
the mixed Hookean-spring theory eq. (8) for x/L — 0,0.47 and 
0.7, from left to right, (b) Corresponding snapshots of the 
filament at points specified in fig. 2 (a). 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

z/L 

(b) 100r 



20 





1 1 1 1 1 r- 




Simuiatior 
eq. m 



_] I I I I I L. 



2 4 6 8 10 12 14 16 

ha/kET 

Fig. 3: (a) Stretching force curves for varying chemical poten- 
tial bias h = 0.0, 7.7, 15.4, obtained at = 0.0124. All other 
parameters are the same as in fig. 2 (a), (b) Threshold force, 
Fc, obtained from fig. 3 (a), plotted as a function of h. Broken 
line is the analytic result, eq. (9), with x ~ 1.32£ as the best 
fit. 



(a). The initial filament shape is coiled, according to the 
experimental conditions [16]. As the extension z increases, 
the helix deforms uniformly like a linear spring. Neglect- 
ing thermal effects, the restoring force -F of a deformed 
helix is on the linear-elasticity level given by [22] 



F = AC 



{kqi cos V' — tqi sin ■0) (Akqi sinip + Ctqi cosip) 

sin-0(ylsin^ t/) -I- Ccos^ V)^ ' 

(7) 

with z/L = cosi/j. For small forces, Hooke's \a.w, 
F ~ Ki,l{z — -LcosV'i), holds, where Km,L = 
dF/dz\,=Lcos^^ - iA7:ye^){Lsm^^Pr„)-\cos''iP^/A + 
sin^ iprn/C)~^ with m — 1,2, is the linear spring con- 
stant [19,22]. The spring constant is inversely propor- 
tional to the arclength of the helix. At a certain critical 
force {F ~ 70 in fig. 2 (a)), the uniform helical shape 
becomes metastable, and a shape transformation occurs 
from the moving end of the filament via nucleation of the 
normal form, see the snapshot C in fig. 2 (b). The fila- 
ment separates into two domains of different pitch angles 
ipi and 1^2, and the tension F suddenly drops down to 
much smaller value; compare B & C in fig. 2. After the 
first transition, subsequent more gradual shape transfor- 
mations at increasing extension z are observed at lower 
forces. Approximating the filament as two serial Hookean 
springs of arclength L — x and x, we obtain a modified 
Hooke's law 



F = 



2.x 



[z — (L — x) cos ipi — X cos 1/^2] ■ (8) 



In fig. 2 (a), eq. (8) is compared to the numerical data 
assuming helix compositions of x/L = 0, 0.47 and 0.7 
(broken lines) and captures reasonably well the response 
of the mixed filament. The main point here is that the 
helical transformation proceeds in subsequent steps. 

A normal helical section with pitch angle -02 appears 
upon stretching via nucleation. Assuming the helix as 
a linear spring, the filament elastic energy (relative to its 
stress-free state) at force F \s E ~ ^F^ / Ki^l. Once a nor- 
mal section of arclength x appears, the tensile force drops 
down to F' at constant extension z. We note that in the 
simulations we almost always observe nucleation accompa- 
nied by two domain walls. The energy of the mixed state 
at the transition thus is E' ~ \ {1/Ki,l-x + l/-fi'2,x)F"^ + 
2hx + 4J. Nucleation proceeds only when E' < E. Us- 
ing F' = (l/i^i,L-. + 1/K2,x)-HF/Ki,l ~ 7), where 
7 = a;(cos'02 — cos^i), and assuming Ki^l-x/ Ki,l — 1 
valid for x/L <C 1, we finally obtain the threshold force Fc 
as a function of x and h: 



Fr. 



-7^2 + ViKi + K2)[K2i^ + 4{hx + 2J)], (9) 



where we have used the shorthand notation Ki^^^^ — Ki 
and K2^x ~ K2- 

We perform stretching simulations for five different val- 
ues of the chemical potential bias, h. All other parameters 
are the same as those in fig. 2 (a). The force-extension 
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Fig. 4: Force-extension curve at velocity V = 0.124, ten times 
faster than fig. 2 (a). Shown is one stretching (red line) contrac- 
tion (green) cycle. The snapshot of the buckling state (left) is 
obtained at stretching z/L — 0.37 upon contraction and yields 
negative forces (highlighted by broken circle). All other pa- 
rameters are the same as used in fig. 2 (a). 

curves for h = 0,7.7,15.4 are shown in fig. 3 (a), from 
which the critical force is obtained as a function of h, 
which is plotted in fig. 3 (b). Since the nucleation do- 
main size X is not determined by our arguments, we fit 
X ~ 19.8a ~ 1.32£ in the comparison of eq. (9) to the 
numerical data in fig. 3 (b), which is also consistent with 
the nucleating domain size extracted from the snapshots B 
and C in fig. 2 (b). The agreement between the numerical 
data and eq. (9) is excellent, verifying the mixed Hookean 
spring scaling model. 

To demonstrate the experimental relevance of our re- 
sults, we put experimental values into eq. (9), namely 
ipi = 76.7°, ^2 = 31.3°, h = 3.4 ^m, £2 = 2.5 ^lul, 
A ~ 3.5 pN /Ltm^, L ~ 7.6 fim [16] and x ~ 1.32^2 ^ 3.35 
/zm from fig. 3 (b). The experimental value of the energy 
difference between the two polymorphs, h, is unknown, 
and we obtain Fc ~ 3 — 5 pN for the range ft, ~ 0.15 — 0.76 
pN ~ 0.04— 0.20A;Br/nm. The experimentally determined 
critical force Fc for the coiled-to-normal transition ranges 
over 3-5 pN [16], quite consistent with our results. The 
direct (and independent) experimental measurement of h 
might be hard, but atomistic simulations [24] may in the 
future check our prediction of h. 

The contraction process involves discontinuous reverse 
polymorphic transformations characterized by sudden in- 
creases of the tensile force; see fig. 2 (b). The reverse 
transitions are in general observed at much lower forces 
and at different extensions, resulting in a large hysteresis. 

Kinetic effects — In fig. 4, we show a force-extension 
curve for the fastest velocity studied by us, V = 0.124, 
with all other parameters as in fig. 2 (a). Notable is the 
negative force in the contraction process at V = 0.124. 
This is due to a buckling of the helix via compression; 
a buckled filament configuration is shown to the left. A 
buckling-induced negative force has also been observed in 
the experiment [16]. 

The domain wall between normal and coiled sections 



moves in response to the pulling of the filament end. We 
plot in fig. 5 (a) the domain wall positions Sd as a function 
of the reduced time Vt/D, where D is the stretching dis- 
tance defined earlier and thus D/V gives the total stretch- 
ing time. The domain boundary position is obtained from 
the discrete cr-profile along the filament. As noted ear- 
lier, we generally observe the occurrence of two domain 
walls. The trajectories in fig. 5 (a) are strongly asymmet- 
ric between the stretch and contraction processes. On the 
stretching side, < Vt/D < 1, trajectories for different 
velocity V collapse. The domain- wall speed dsd/dt upon 
stretching approximately follows from a simple geometric 
argument as 



dt cos 1/12 — cos -01 

To see this, note that the lengths of the two helices are Sd 
and L—Sd- Neglecting helix deformations from its intrinsic 
shape, the end-point position follows as 2; ~ SdCOSil)i -\- 
{L — Sd) cos V'2: which upon inversion yields Eq.(lO) where 
V = dz /dt is the imposed stretching speed. In the reduced 
units in fig. 5 (a), this gives the fine of slope (cosf/'2 — 
cosipi)^^ independent of V, which is shown as a broken 
line, in reasonable agreement with the numerical data at 
early times. Deviations at later stages are connected to 
nonlinear filament deformations. 

To quantify rate-dependent dissipation effects, we cal- 
culate the work done between zi = 0.3L and Z2 = 0.8L, 
AW — J^^ dzF, in both the stretching and contraction 
processes, which is shown in fig. 5 (b) as a function of V. 
The difference between the two works, which is a mea- 
sure of dissipation, becomes larger as V increases. Es- 
pecially for small V, AW is strongly influenced by the 
stretching rate V, consistent with the experimental ob- 
servations [16]. One trivial contribution to dissipation 
is due to solvent friction. Assuming the fllament to be 
homogeneously stretched and on the free-draining level 
(as appropriate for our simulations), the friction work is 
AWf — Trr]aN{z2 — zi)V, which is shown in flg. 5 (b) as 
broken lines. As one can see, the solvent friction does not 
account for the entire dissipation, the difference is due to 
elastic dissipation within the filament. In the quasi-static 
limit y — > 0, AW for stretching and contraction should 
become equal. Although a direct numerical check of this 
is difficult because of prohibitively long simulation times, 
a simple scaling estimate is helpful: Assuming that the 
helix takes a perfect normal form (all a = —1) at high 
stretching z = 0.8L without any elastic deformation, and 
a perfect coiled form with all ti = 1 at 2: = 0.3L, the total 
elastic energy change is simply AE = 2hL which deter- 
mines the work AW in the quasistatic limit at = 0. In 
rescaled units, this leads to AW/ksT ~ 2hN = 924 for 
h = 7.7, which in fig. 5 (b) is indicated with a solid hori- 
zontal line and is flanked by the stretching and contraction 
data as expected. In the large V region, i.e. for V > 0.06 
in flg. 5 (b), which probably is beyond experimentally ac- 
cessible speeds, the force response seems less affected by 
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Fig. 5: (a) Trajectories of the rescaled domain-wall positions 
Sd/L as a function of the reduced time Vt/D for varying ve- 
locities V . The slope of the broken line is (cos 1/^2 — cos'i/'i)"^ 
from eq. (10). (b) Mechanical work done between z\ = 0.3L 
and 22 = 0.8L, AVK/fcsT, as a function of V in the stretching 
(square) and contraction (circle) processes. Each data point is 
the average over three statistically independent runs. (Error 
bars are smaller than the symbols.) The solid line indicates the 
equilibrium work, and the broken lines are the friction work for 
stretching and contraction processes (see text) . Parameters are 
the same as in fig. 2 (a), unless otherwise stated. 



the stretching rate V and the apparent rate-dependence is 
largely due to solvent friction. 

Summary — In this letter, wc have introduced a 
bistable helical filament model that accounts for differ- 
ent elastic monomeric states by a discrete Ising-like spin 
variable along the arclcngth. Employing hybrid Brownian 
dynamics Monte-Carlo simulations, we have studied the 
kinetic effects of uniaxial stretching on the helical shape. 
The main results arc: (i) The helical filament imdcrgoes 
reversible polymorphic transformations in response to an 
externally applied uniaxial tension, consistent with the 
optical-tweezer experiments of Salmonella fiagcUar fila- 
ments, (ii) Experimental force-extension curves includ- 
ing sawtooth pattern and large hysteresis are qualitatively 
reproduced, (iii) A simple analytic argument based on 
a mixed Hookean spring model for the critical force of 
the polymorph induction provides good agreement with 
numerical and experimental data, (iv) Force-extension 
curves are strongly pulling rate-dependent, in particular 
a buckling-induced negative force is observed, all consis- 
tent with the experimental observations. 



Generalization of our model to include all possible poly- 
morphic flagellar states is straightforward but includes 
many at present unknown parameters. Applying our 
model to torque-induced chirality transformations of flag- 
ellar filaments, directly relevant to bacterial motility in 
vivo, is in progress [25]. 
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